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ABSTRACT 


rLeceumoepmbocs caused by lrradiation of a model gallium 
arsenide solar cell with a Carbon Dioxide TEA Laser 
were examined using power densities of the order of 
iho” watts/cm’. The pressure pulses were monitored with a 
Sandia type quartz pressure gauge. It was discovered 
that the relatively low power densities used were capable 
of removing the silicon dioxide antireflective and gold 
contact layers of the solar cell after only a few shots 
of the laser. An exponential relationship between the 
initial thicknesses of these layers and the pressure 
pulse generated in the gallium arsenide substrate was 
indicated for gold contact layers of less than 5000 A 
thickness. Evidence was found that the principal pressure 
generation Mechanism is thermomechanical. Gold films of 
thickness greater than 5000 A were found to be able to 


absorb the power densities used with no apparent damage. 
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I. INTRODUCTION 


The development of solar cells as versatile and dependa- 
ble sources of power for satellites is well-known. Virtually 
every long-range space probe launched by the United States 
or the Soviet Union in the last decade has been powered by 
photovoltaic devices. Solar cells provided the power that 
Operated the television cameras of the Ranger, Surveyor and 
Lunar Orbiter spacecraft that made possible detailed photo- 
meaphic Studies Of the surface of the moon. These devices 
provided the power for the Mariner probes of the planet 
Venus. The recent Skylab vehicles also utilized solar cells 
in their electrical generation systems. The application of 
solar cell technology has led to the development of the 
systems of communications, navigational and meteorological 
satellites currently being employed by the United States 
and other countries [Ref. 1]. 

The general construction and characteristics of solar 
cells are treated in detail by several authors. A good 
review of these topics is contained in Direct Energy 
Conversion by S. W. Angrist [Ref. 2]. Of particular note in 
this treatment is the requirement that for maximum theoretical 
conversion efficiency the solar cell material should have an 
energy gap of between 1.1 and 2.3 electron volts. Gallium 
Arsenide (GaAs) with an energy gap of 1.35 ev. satisfies 
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A significant amount of research has been undertaken 
previously to ascertain the effects of laser irradiation 
on solar cells and phototransistors. A study of irradiation 
of phosphorus-doped silicon solar cells by a pulsed laser 
performed by Giuliani [Ref. 3] identified two permanent 
damage effects: Below laser energy densities of 5 joules/cm” 
no permanent damage to the phototransistor was observed. 
Irradiation at these low energy densities was found to 
increase the photoconductive response of the device, however. 
At energy densities above 10 joules/cm< permanent changes in 
the electrical characteristics in the form of increased 
dark current were observed. Examination of phototransistors 
damaged by laser energy densities in this range enabled 
Marquardt, Giuliani and Fraser [Ref. 4] to explain this 
increased dark current. Evidence was found that the melting 
point of the silicon had been exceeded. A model for phos- 
phorus ion migration was developed by treating the phosphorus 
as a dilute impurity in a melted silicon matrix and then 
analyzing the migration of the phosphorus ions as a diffusion 
problem. The subsequent refreezing of the silicon then traps 
the ions in their new distribution and causes the increased 
dark current observed. 

As evidenced by the above work a mechanism explaining 
the changes in the electrical properties of phototransistors 
under laser irradiation has been developed. However, the 
mechanical damage done to solar cells by laser radiation 


has not been so thoroughly examined. By investigation of 
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the momentum delivered to a GaAs solar cell by laser irradia- 
tion it was the intent of the current investigation to amplify 
the knowledge in this area. A model of a GaAs solar cell 

was constructed as shown in Figure 1. The thickness of the 
gold film was varied and a measurement of the transmittance 
'of the gold film was made at the wavelength of the Carbon 
: Dioxide laser used (10.6 um). This was done to determine 
whether the gold coating would absorb a significant amount 
of the incident laser energy. The solar cell model was 

then constructed and irradiated with light from a co, TEA 
laser. The impulse transmitted to the GaAs was measured as 

a function of time and of the thicknesses of the gold and 
S10, film using a Sandia-type quartz gauge [Refs. 5 and 6]. 
Finally, an experiment was performed to determine the energy 


threshold for removal of the quartz film by the incident 


laser beam. 
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II, PREVIOUS WORK 


The absorption of laser radiation by a great variety of 
materials under a large range of conditions has been examined 
by several experimenters. The type and mechanism of damage 
is dependent upon whether the target material is opaque or 
transparent to the incident laser light. 

Damage to transparent materials by laser radiation has 
been the subject of much experimental and theoretical work. 
Observations of the damage itself are plentiful but a clear 
understanding of the mechanism of damage has yet to be 
developed. 

The main mechanisms suggested for the laser-induced 
damage to transparent materials are [Ref. 7]: 

a. Production of phonons (hypersound) in a stimulated 
Brillouin scattering process with the phonons being produced 
in an intensity sufficient to lead to a fracture of the 
material. 

b. Absorption of light by defects in the material with 
resultant thermal shock and microplasma production. Using 
a ruby laser with 20 psec to 20 nsec pulse durations and 
depositing laser energy densities of from 3.5 joules/cm? to 
IL ea, joules/cm” on the target, Milam, et al. [Ref. 8] demon- 
Strated that heating of inclusions ranging in size from 0.1 


to 5.0 um is the limiting mechanism for damage in TiO, and 


2 
Zr0., films deposited on quartz substrates. A relationship 
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between the length of the laser pulse and the size of the 
inclusion heated was also developed. Short laser pulses 
can more easily heat small inclusions and long pulses 
preferentially heat large inclusions. It was shown by 
Hopper and Uhlmann [Ref. 9] that, in the first approxima- 
tion, the inclusion that is most readily heated by a laser 
pulse of duration T) has a radius proportional to © et 

c. The electron avalanche phenomenon. This phenomenon 
occurs at the surface of a transparent dielectric [Ref. 10]. 
"Free'’ electrons at the surface of the dielectric are 
accelerated by the electric field of the laser beam sufficient- 
ly to produce ionizing collisions with the atoms in the 
surface layer thereby producing more free electrons. When 
the lattice of the solid acquires energy from the free 
electrons so produced its temperature rises. The surface 
matter then melts, boils and evaporates, forming a partially 
10nized high-density vapor layer above the surface. This 
vapor layer absorbs laser light strongly and the absorbed 
light energy continues the ionization process. The mechani- 
cal damage resultant at the surface is primarily due to the 
plasma formed as a result of an optical discharge (thermal 
ionization of all the evaporated matter) and a shock wave 
generated by this discharge. 

Damage produced in transparent materials by high-energy 
laser beams can range over a wide field of effects extending 


from microcracks or melted voids in the interior of the 
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material to craters and cracks on the surface of the material 
[Ref. 11]. 

The damage effects in transparent materials are highly 
power-dependent. If the incident laser radiation is below 
a certain power threshold, the laser beam passes through the 
material with no effect. However, if this threshold is 
exceeded, the laser beam is absorbed and damage of one form 
Or another takes place [Refs. 12, 13 and 14]. For sodium 
chloride crystals irradiated with a Q-switched ruby laser 
this threshold occurs at 2.0 GW/cm?: fot oer re om, 1 oxemcde 
glass the threshold is at 470 GW/em? [Ref. 12]. Wang [Ref. 14] 
has measured the damage threshold for copper reflectors 


irradiated with 10.6 um radiation from a CO. TEA laser to be 


2 
150 MW/cm?. Experiments [Ref. 11] have shown that the 
threshold for internal damage occurs at lower levels of 
irradiation than surface damage. 

The absorption of laser emsiensy by substances opaque to 
the wavelength of the laser used (absorbers) has also been 
the subject of considerable investigation. The variance of 
energy absorption with irradiance and duration of exposure 
shows three distinct regimes differing in the phenomena 
induced in the material by the laser beam. 

The first of these regions is heating of the material 
Without phase change. This results in a temperature rise 
in the material. This temperature rise may lead to simple 
thermal expansion of the material if the deposition of energy 


1s slow enough. However, if the laser energy is deposited 
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rapidly enough, say in a time t comparable to or less than 
the time required for an elastic wave of velocity v to 
cross the area irradiated, a pressure wave can be generated. 


Phat is, if 


t < a/v (1) 


where a is the radius of the area irradiated. In this case 
the irradiation takes place so quickly that the material 
remains "inertially clamped" during the absorption process. 
Hence, the volume of the irradiated portion of the material 
reMains constant over the time immediately following the 
laser irradiation. A laser pulse of duration aaa sec. 
Saedoitles this condition [Ref. 15]. 

The next region considered is laser-induced melting of 
the target without vaporization. In this region the damage 
caused in the target consists of melting of target material 
to a depth dependent upon the power density of the laser 
beam irradiating the target and upon physical properties, 
such as the heat of fusion, of the target material. The 
optimum pulse length required for this effect ranges from 
50 usec to 5 msec and is much longer than the pulse lengths 
in the nanosecond range that were produced with the lasers 
used in the current investigation [Ref. 16]. 

The final region in which absorption by opaque materials 
takes place in laser-induced vaporization. For the purposes 


of the experiment performed this was the most significant 
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region. The strong energy fluxes concentrated over the small 
area of the focal spot of a laser beam produce such intense 
evaporation that significant pressures can be generated in 
the absorbing material by the recoil of matter ejected by 
evaporation [Refs. 17 and 18]. When the irradiation takes 
place in a vacuum a hot, dense plasma is formed above the 
material and the recoil of the particles in the plasma 
transfers momentum to the target. In an experiment in which 
aluminum targets were irradiated with a Neodymium-glass 
laser having an output between 4 and 11.5 joules, Krehl, 
Cooper and Schwirzke [Ref. 19] demonstrated that of the total 
mass ejected by the laser beam from the target about ten 
per cent was ionized and that this small ionized component 
was responsible for almost 80% of the momentum transferred 
to the target. A simple calculation by Askar'yan and 
Moroz [Ref. 17] shows that the recoil pressure of the ejected 
matter 1S many times the radiation pressure of the laser beam. 
Pirri [Ref. 20] has developed a comprehensive theory for 
the momentum transfer to a surface from a laser-supported 
detonation wave initiated in the air above the target surface 
by a laser beam. Although this theory is not strictly 
applicable to the vacuum portion of the current investigation, 
it is described here.in some detail as an excellent descrip- 
tion of the mechanism of momentum transfer to a surface by 
a laser beam. The model is based upon the physical processes 


which occur when a high-power axisymmetric laser beam 
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interacts with a surface in the presence of an atmosphere. 
Figure 2 illustrates this model. 

am Ty 1s the laser pulse length, at some time t << The 
on the order of 50 nsec, a laser-supported detonation wave 
is initiated above the target surface by breakdown of the 
small amount of surface material vaporized during this time. 
Once the detonation wave forms, the vaporization stops and 
the force delivered subsequently is a result of shocking the 
air above the surface to a high pressure. 

This high pressure in the laser-shocked air means that 
there 1S a significant impulse delivered to the target 
after the laser pulse ends as the pressure relaxes to its 
ambient value. Hence, the dominant time scale for the 
momentum transfer process is Bas the time over which the 
impulse is delivered to the target surface rather than just 
the laser pulse length Tp: 
Initially, the flow of the hot, dense plasma from the 


surface is one-dimensional. At some later time, the 


"2D? 
flow becomes two-dimensional due to the thermal expansion 
of the plasma. This two-dimensional time is of the order 
of the time it takes the leading edge of the plasma to 


travel a laser spot diameter [Ref. 21]. This can be defined 


as 


S 
ee g dr/V (2) 
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where D. is the diameter of the laser spot and V is the 
velocity of the leading edge of the plasma. 

For the entire pressure-time history of the momentum 
transfer to be described by one-dimensional analysis the 
total time over which the impulse is delivered must be less 


than the two-dimensional time t... That is T, Saas 


2D Ze 


In one dimension the detonation moves away from the 
surface with a velocity Vow? given by 


i ae 


Vow = [2(y*-1)E,/9,] (3) 


where y is the ratio of specific heats of the air behind 
‘the detonation wave, E, is the laser flux and 9, is the 
ambient density in front of the wave. The pressure behind 


the Wave is 


Pow = PoVpw /Cyth) - (4) 


The surface pressure on the surface while the laser is on 


and the flow is one-dimensional is 


Psrp 7 (Cyst) /2y]2 Ye (5) 


The impulse per unit area at a point on the surface, l, 


is the sum of the contributions to the impulse transferred 
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while the laser is on and after the laser pulse has terminated. 
It is shown in later sections of this report that a one- 
dimensional analysis is justified for the conditions of the 


current investigation so only the case where t Will 


ye OOD 
be considered. 


While the laser pulse is on, the impulse transferred to 


the surface is 


After the laser pulse is off the shock expansion flow field 
relaxes as a planar blast wave until Top is reached. The 


pressure is then 
eee) = Beet sry! (tT. < t < Toy) (7) 
S JU as - ~— “2D 


The impulse delivered to the surface is then 


Cc 
a5 
= p lg Ee 
een Sop. te yf dt/t 
0 


= a ie 
= 2s Be {1 3[(sip/Pscty? - lj} (0 < t < Typ) (8) 


The total impulse is now obtained by setting p(t) equal to 


ee the ambient pressure in front of the detonation wave. 


To 





Then, by substituting equations 3 through 5 into equation 8 


and by assuming Perna. oo? it can be shown that 


for tT, S Typ, Where V, is the ambient speed of sound. The 
impulse delivered to the surface is, therefore, proportional 
to the speed of sound in the ambient gas. 

Experiments performed by Pirri, Schier and Northam [Ref. 
22] showed the theory developed above to predict satis- 
factorily the momentum transfer to a surface irradiated with 
a high-power pulsed CO. laser in an average flux regime 
between 4 and 200 MW/cm?. 

Although Pirri's theory presents a comprehensive treatment 
of momentum transfer to a surface by a laser it is questionable 
whether it is completely applicable to the conditions used in 
BitoeCcULbrenc Investigation. According to Harrison [Ref. 23], 
in order for a laser-supported detonation wave to be formed 
laser power densities of the order of ine to une watts/cm- 
are required in an atmosphere of pressures above 76 torr. 
Portions of the present investigation were carried out at 
pressures of the order of ae COG tecnese low pressures 
the impulse delivered to the target as a result of shocking 


the air above the target is not present. If this contribution 


to the total impulse is neglected, the impulse delivered to 
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the target is 


I = Point (t < tJ (10) 
This 1S just the expression for impulse per unit area 
delivered to a target in vacuum reported by Gregg and 
Thomas [Ref. 24]. 

A more suitable equation for predicting the maximum 
pressure developed in a target is given by Anderholm 
(Ref. 25]. Here the target is treated as a perfect gas 
bounded by material with acoustic impedance Z. The maximum 
pressure developed is 


1/2 (11) 


PAX = (ZE,/3t,,) 
where ES is the laser energy density as before and Ty is the 


laser pulse width. The acoustic impedance is given by 
(12) 


with p being the mass density and V, the bulk sound velocity 
in the target. Verification of the maximum pressure 
predicted by this equation was obtained using a ruby laser 

to generate a pressure of 34 kbar in a quartz target. The 
value for maximum pressure of 29 kbar predicted by equation ll 


agreed satisfactorily with the experimentally measured value. 
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Irradiation of opaque targets can generate very high 
pressures in the target. Several experiments have been 
performed by various authors to verify this fact. By 
irradiating a 0.66 mm thick aluminum target with a pulse of 


power density 2000 watts/cm- from a CO, laser, McMordie 


2 
and Roberts [Ref. 26] measured pulses of up to 2 Bar. As 
noted above pressures in the kbar range were measured by 
Anderholm. By coating 1 mm thick targets with a water and 
paint coating and then irradiating the targets with light 
from a Q-switched Neodymium:glass laser with a laser 
fluence of 25 joules/cm“, Stresses in the kilobar ranges 
were measured by Fox [Ref. 27]. 

If the absorbing material is confined in some way, as 
in the experiment by Fox described above, the mechanism 
generating the stresses in the target is somewhat changed. 

Confinement of the absorbing material between two trans- 
parent films causes a confined high-temperature plasma to 
be produced when the target is irradiated by a laser. This 
confined plasma drives a shock wave into the transparent 
media [Ref. 28]. The amplitude of this shock can be of the 
order of tens of kilobars. 

Covering the absorbing material with a thin coating 
and then illuminating it with a high-power laser brings a 
different effect into play. The thin film is partially 
vaporized as a result of heating by thermal conduction from 


the laser-heated absorber surface. The resulting expansion 
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of gas and partial confinement at the absorber/thin film 
interface causes stress waves to be generated in the 
absorbing material [Ref. 29]. Experiments performed by 
Magee, Armistead and Krehl [Ref. 29] and by Fox [Ref. 27] 
have shown that the volatility of the coating has a signifi- 
cant effect upon the amplitude of the stress wave produced. 
Magee noted that as the curing time of the silver paint 
used as a coating for his targets was increased the magni- 
tude of the stress waves produced in the target was reduced. 
This result was attributed to a change in the transmittance 
of the paint as its butyl acetate binder evaporated. In 

the experiment by Fox reported above [Ref. 26] it was 
observed that the coating of the target with a mere drop of 
water Significantly increased the amplitude of the stress 
wave produced. For an aluminum target 0.17 mm thick it 

was noted that backface spallation occurred at laser fluences 
around 310 joules/cm- when the target was water-coated. 

For an uncoated target fluences of up to 500 joules/cm” 
produced no spallation. If the target was also coated with 
black paint, the threshold for backface spallation was 


further reduced to 200 joules/cm*. 


Zo 





IIIT. EXPERIMENTAL DESIGN 


A. EQUIPMENT 

A Carbon Dioxide (CO, ) Transverse Excitation at 
Atmospheric Pressure (TEA) laser was utilized throughout this 
experiment. A description of the laser is given by Strickland 
[Ref. 30], Behrens [Ref. 31] and Stevenson [Ref. 32]. The 
laser output was reflected off a 1m focal length copper 
mirror through a germanium lens onto the solar cell model. 
The solar cell was, in turn, mounted in a test chamber that 
could be maintained either at atmospheric pressure or at a 


waeuem of 10°” 


torr. depending upon the conditions desired. 
The solar cell model was mounted on a Sandia-type Quartz 
Pressure Gauge which enabled the pressure at the back face 
of the model to be monitored on a Tektronix 7704 oscillo- 
scope. The energy output of the laser was monitored using 
a Ballistic Thermopile. A focal spot size of approximately 
1.0 cm was maintained on the solar cell model by means of the 
germanium lens. A schematic of this experimental arrangement 
1s shown in Figure 3. 

she co, TEA Laser 

The Naval Postgraduate School's CO, TEA Laser System 

consists of a resonant cavity, a Tachisto Marx-type generator, 


a pasSive electrical component pulse-shaping network, a direct 


current power source and a gas flow system. 
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The resonant cavity consists of a gas-filled lucite 
Pere weconstmucted of 2.54 cm thick lucite, 101.6 cm long, 
30.5 cm wide and 15.3 cm high, and two mirrors. At each 
end of the. box a Sodium Chloride (NaCl) crystal is mounted 
to serve as a Brewster window. The back end of the resonant 
cavity is a 10.2 cm diameter polished copper mirror having 
pom Goecal length. The front minrorsis an 11.4 cm 
diameter plane parallel germanium mirror which transmits 
approximately 47% of the radiation striking it from inside 
the cavity, thus providing an outlet for the laser radiation. 

Inside the lucite box is the electrode assembly which 
consists of a 2.0 cm thick aluminum plate located at the 
bottom of the box and a cathode assembly at the top. The 
plate is cut to fit the inside of the box and serves as 
the anode of the system. On the inside top of the box is 
the cathode. It is composed of 131 aluminum blades separated 
by lucite spacers and electrically connected by two threaded 
brass rods. The trigger assembly consists of 130 glass- 
encased nichrome wires set between the aluminum cathode 
plates with each wire set at the same height as the edge of 
the cathode blade. A gap of 5.8 cm is maintained between 
the anode wires and the anode plate [Ref. 32]. Some 
characteristics of the laser are given in Table I of 
Appendix A. 

Selective excitation of the cO., upper laser level 


(the 001 level) in a TEA laser is collision dominated. An 
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electric discharge across the gap between the anode and 
cathode injects energy into the laser cavity which creates 
a uniform flux of electrons having an average energy of about 
one electron volt. A large number of inelastic collisions 
then take place between the CO. molecule in the ground 
State and the excited electrons. This causes the carbon 
dioxide molecule to be excited to various vibrational levels. 
Resonances exist in the carbon dioxide excitation cross- 
section for electron energies of 0.3, 0.6, and 0.9 electron 
volts. These resonances coincide with the vibrational 
ladder of energy levels of the carbon dioxide molecule and 
hence preferential excitation to these levels occurs. Thus 
an electron energy in the glow discharge near one electron 
volt excites the desired energy level ladder. Collisions 
between excited and unexcited carbon dioxide molecules thus 
populate the upper laser level very efficiently. The high 
efficiency arises because the energy levels are equally 
spaced and little energy is lost during transitions up 
or down the vibrational ladder. This high excitation effi- 
ciency, coupled with a quantum efficiency of about 41%, 
enables an operational efficiency of around 25% to be 
obtained. (Quantum efficiency is defined as the ratio of 
the laser photon energy to the upper level excitation 
energy.) 

The presence of nitrogen (N.) in the laser atmosphere 


provides a complementary excitation process. Nitrogen 
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molecules are also excited to a higher vibrational energy 
level by electron impact. The lowest excited vibrational 
energy levels of nitrogen and the upper laser level of the 
carbon dioxide molecule are closely matched so an efficient 
transfer of energy occurs through collision of excited 
nitrogen molecules with ground state carbon dioxide molecules. 

Once a carbon dioxide molecule is selectively 
excited to the upper laser (001) level, it can emit a 10.6 um 
photon as it decays, or is stimulated, to the lower laser 
level (010). It must be quickly removed from this level to 
maintain the population inversion necessary for lasing. This 
1s accomplished by adding a light gas, such as Helium (He), 
with a collision rate much larger than carbon dioxide to the 
laser atmosphere to increase the rate of depopulation of the 
lower laser level to the ground state. 

The NPS CO. laser is operated with an atmosphere of 
carbon dioxide, nitrogen and helium in a mixture of 90/5/5 
per cent and at a flow rate of about one-half liter per 
second. 

Conditions for lasing are induced by electrical 
breakdown between the cathode and anode located in the laser 
box. This produces the direct electron collision and 
molecule-molecule energy transfer process that results in 
the carbon dioxide population inversion. The laser is 
operated at atmospheric pressure rather than at low pressure 
in order to make use of the high molecular density and 


consequently greater energy available. 
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The electrical power required for excitation of the 
laser is produced by a three-stage Marx generator supplied 
with a positive power voltage of up to 22 kilovolts and 15 
milliamperes by an unregulated d.c. power supply. Due to the 
voltage multiplication of the Marx generator voltages of up 
to 66 kilovolts can be generated across the laser electrodes. 
Pulse-shaping 1s required to obtain the uniform glow discharge 
needed for gas pumping. This is accomplished through the use 
of a network of passive circuit elements that controls the 
voltage rise time and decay time constants of the input 
energy and the Marx generator. This network is discussed 
moe net. 31. 

2. Solar Cell Model 

Reference 33 describes the fabrication of a Gallium 

Arsenide Solar Cell. The Gallium Arsenide (GaAs) substrate 
is coated with a continuous platinum (Pt) film approximately 
40 Angstroms thick to establish the barrier contact height. 
This platinum film is then overlaid with a gold (Au) contact 
several thousands of Angstroms thick. Since a gold layer 
of this thickness is essentially opaque to visible light, 
a large area solar cell such as is used on satellites, 
utilizes a grid structure for the gold contact layer. Such 
a grid structure is illustrated in Figure 4. The entire 
solar cell is then coated with an antireflective coating of 
Silicon dioxide (Si0,). 

The model of the gallium arsenide solar cell con- 


structed for the current investigation was designed to 
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approximate the cross-Section of the cell at one of the 
grid lines. An illustration of this model is given in 
Figure 1. The model consisted of a gallium arsenide sub- 
Strate approximately 3.0 mm thick and 1.0 cm in diameter. 
A gold film of thickness varying up to 5000 Angstroms was 
vacuum deposited on this substrate. <A silicon dioxide film 
of thickness up to 3000 Angstroms was then deposited onto 
the gold film. The model was then mounted on a Sandia quartz 
gauge uSing a thin layer of epoxy cement. 

Since the gold layer is thick enough so that all 
of the laser radiation incident upon it is absorbed, the 
momentum coupling between the gold and the platinum layers 
would be entirely mechanical. It [Ref. 34] has been stated 
that a film of thickness on the order of 50 Angstroms is 
invisible as far as mechanical coupling is concerned. 
Therefore the thin platinum film was not included on the 
solar cell model used in the current investigation. 

3. Sandia Quartz Gauge 

The piezoelectric quartz gauges used in this experi- 
ment were provided by Stanford Research Institute and were 
of the type developed by Graham [Ref. 5] at Sandia Labora- 
tories. This type of pressure gauge has a resolution on the 
order of a few nanoseconds. 

The operation of the gauge is based on the fact 
that when a force is applied on a piezoelectric crystal, 


such as quartz, it results in the formation of a stress 
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fmrmin the crystal that distorts the ionic structure of the 
crystal. This distortion results in an increase in the 
polarization of the sample. 

It 1s well-known from electromagnetic theory that 


the electric displacement vector is given by 


= > 
De=-ck + P (13) 
san ; > 
Mionewoeis the plezoelectric polarization, E is the electric 
field and € is the dielectric permittivity. If a one- 
dimensional strain 1S assumed, the piezoelectric equation 


becomes 


= Nal toe (14) 
where T is the stress and do4 1s the piezoelectric stress 
constant. If it is assumed that the piezoelectric polariza- 
tion is directly proportional to the stress with a coefficient 
f that is independent of time and stress over a given stress 


range, ee becomes 


P_ = fo, (x) (15) 


where og, (x) is the x-component of the stress. For an X-cut 


quartz crystal [Ref. 5] 
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ad. = 2.25° 10°14 


21 coulomb/newton 


| ee 06 108 5 


farad/meter. 
For the static case, the open circuit voltage is given by 


[Ref. 35] 


g = d,,/e = 0.55 2£F/Tw (volt-meter/newton) (16) 


and the short-circuited charge for a given applied force is 


Q= ds4 Bet er, 
Where F is the applied force, considered positive for an 
extensional force, 2 1s the crystal length, w is the crystal 
width and t is the crystal thickness. 

The development given above is adequate for low- 
frequency impulses. However, in the current investigation 
the analysis is complicated because the laser pulse produces 
a stress wave with a short wavelength. 

Graham, et al. [Ref. 5] have analyzed the short- 
circuit piezoelectric current induced by a rapidly changing 
impulsive load applied to the quartz. This load propagates 
as a stress wave along the x-axis of a quartz disc. The 
displacement current generated by this stress wave is 


i SA dD/dt (18) 
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mere A 1S the electrical area of the disc and D is the 
displacement vector. 
From equation 14 and assuming one-dimensional stress 


in the eect e ect On) 1t is seen that 


n n q 
ne D(x) dx = Fi P(x) dx +e f EB (x) dx (19) 


@) O O 


where 2 1S the thickness of the quartz diSc. 
Assuming a short circuit between the electrodes 
and that the crystal has a constant permittivity, the last 


integral in equation (19) becomes 


Q 
f cE (x) dx = 0 (20) 


O 


The conductivity of the quartz is essentially zero, 


SO 


Ye 


1/2 f P(x) dx =D, . (21) 


O 
From equation 15 it is known that 


Ee = jee) (oe) 
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Pe it 1S assumed that the quartz is linearly elastic it can 


be stated that 
Onexe t) = o(x-u,t) 22 } 


where u. 1s the wave propagation velocity. Therefore, from 


equation 18, we have 


Q 
A d/dt [1/2 I fo(x-u.t) dx] 
O 


a0 


x 
- -fu A/2 fo aniox)/ oxwdx 
O 
=f u. A Cai - O)/k ESD 


where o 5 1s the x-component stress at the input electrode 
and o. is the x-component at the rear electrode. 
For times less than the time required for the stress 


wave to travel the length of the quartz, a = 0 and 


ln = f A u. oof . (24) 
Therefore, for times less than the wave transit time, the 
piezoelectric current generated in the quartz crystal is 
directly proportional to the stress generated in the quartz 


by the laser pulse. The gauge used in this investigation 
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has a time resolution on the order of a few nanoseconds and, 
using a 50 ohm shunt resistor, has a sensitivity of 1.09 
mV/bar. The gauge is illustrated in Figure 5. 

This gauge employed a guard ring configuration as 
described by Jones [Ref. 6] and illustrated in Figure 6. 

As is shown, the inner portion of the X-cut quartz disc is 
electrically isolated from the outer portion by separating 

the vapor-coated electrode into two regions. This guard 

ring configuration eliminates field distortions. At the outer 
edge of the disc, the discontinuity in electric potential 

and dielectric permitivity of the AG aig produces electric 
field fringing analogous to that found in parallel-plate 
capacitors. The guard ring restricts observation to the 
central region of the quartz disc where there is negligible 
Pane ine . 

This guard ring also aids in the elimination of 
distortions caused by unloading waves. As the stress wave 
propagates through the disc in the axial direction, boundary 
conditions generate shear and dilational waves immediately 
bolmendetheswayvetront at the lateral edge of the disc. These 
unloading waves then propagate laterally inward from the 
edge and the crystal is no longer in a state of one-dimensional 
Strain. The central region of the disc will be one-dimensional 
only for the first wave transit time if the width of the 
center electrode is such that the unloading wave does not 


reach the central region during the first wave transit time. 
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For proper monitoring of the pressure pulse, the 
propagation time of the pressure pulse in the quartz must 
be larger than the half-width of the laser pulse. For the 
250 ns pulse width of the co, TEA laser and a speed of sound 
in quartz of 5700 meters per second, this requires that the 
gauges used in the cO., laser must have quartz discs at least 
1.44 mm thick. The gauges used had quartz discs 3.0 mm 
thick so this requirement was satisfied. 

Care also had to be taken that a pulse reflected at 
the rear of the gallium arsenide target sample did not reach 
the gauge during the loading phase. This required that the 
target sample thickness, d_, had to satisfy [Ref. 29] 


d. >t 


s eae (25) 


q 
where ty is the propagation time of a pressure (sound) pulse in 
the quartz plate and We is the speed of sound in the gallium 
arsenide. For a sound speed of 5000 meters per second in 
gallium arsenide and the 500 ns transit time in the quartz 
disc calculated from the data for quartz given above, this 
required that d. be greater than 1.25 mm. The 3.0 mm thickness 
of the gallium arsenide samples used satisfied this requirement. 
A final requirement for proper operation of the quartz 
gauge was that the stress be essentially uniaxial. To ensure 
this a laser spot diameter of slightly greater than the 1.0 
cm diameter of the gallium arsenide sample was used 


throughout the experiment. 
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4. Ballistic Thermopile 


The total energy output of the CO, TEA laser was 


2 
monitored by directing the laser beam into a ballistic 
thermopile. The ballistic thermopile is a calorimetric 

device that measures the temperature rise due to absorbed 
radiation. It is composed of two nickel-plated cones 
(receiver and reference) and two series-connected iron- 
constantan thermocouples. These thermocouples are situated 

so that the hot junctions are attached to the receiver cone 
and the cold junctions to the reference cone. Radiant 

energy is directed into the receiver cone where it is almost 
totally absorbed. The temperature difference between the 
receiver cone and reference cone generates an emf, which is 
monitored by an externally connected microvoltmeter. The 

peak meter reading 1s linearly proportional to the total input 
energy with a constant of proportionality equal to 44.9 uV/ 
joule. Due to the high reflectivity in the cone, the pol- 
ished surface resists the destructive effects of the high 

peak power pulses found in laser plasma research. Some 
characteristics of the HADRON Ballistic Thermopile used in 
this investigation are given in Table II of Appendix A. 


5. Test Chamber for CO, Laser System 


The test chamber utilized was constructed of a 
rectangular aluminum block 11.4 cm long, 8.9 cm wide and 
10.0 cm high through which a 6.4 cm diameter hole had been 


drilled from the top to the bottom. Five access ports were 
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provided, a 4.0 cm diameter port in each side and a 6.4 cm 
diameter port on the top of the chamber. An overhead view 
of the test chamber is shown in Figure 7. 

Laser light was introduced into the chamber through 
Port A. Ports B and C were covered with plexiglass windows 
to enable the system to be placed under vacuum by means of 
a vacuum system connected at Port D. The test chamber 
was mounted on an aluminum bed plate and that was adjustable 
in height by means of bolts mounted at each corner of the 


pace’. 
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1. Determination of Peak Pressure in the Gallium 
Arsenide Substrate 


Equation 10 provided a method whereby the maximum 
pressure developed in a transparent material by laser 
irradiation can be calculated. However, the problem 
presented by the generation of a pressure pulse in the GaAs 
layer of the solar cell model used in the current investi- 
gation is complicated by the fact that, although GaAs is 
essentially transparent to the 10.6 um radiation used, the 
gold and silicon dioxide coatings are strongly absorbing 
at this wavelength. Also, the thickness of these two films 
was such that the incident laser radiation was totally 
absorbed prior to reaching the GaAs substrate. Hence the 
problem of predicting the magnitude of the pressure pulse 


in the GaAs was reduced to that of predicting the pressure 
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in the gold-silicon dioxide combination layer and then 
computing a correction factor for the mechanical coupling 
of the pressure pulse across the gold-gallium arsenide 
interface. 

By consideration of the conservation of momentum 
across the interface between two materials of different 
acoustical properties Krehl, et al. [Ref. 19] developed an 
equation for this correction factor. Employing this result, 
the magnitude of the pressure in the gallium arsenide 
transmitted by the gold layer is 

P 


[Oe (26) 


Au Fr Pcaas “Gass °GaAs id We AnD GaAs*GaAs 


where op is the initial density and v is the elastic wave 
velocity in the medium. Using the values of o and v for 
gallium arsenide and gold given in Table III of Appendix A 


we obtain 


Deane’ | At = Use (27) 
A similar correction must be made to the pressure measured 
by the quartz gauge to obtain the pressure in the gallium 
arsenide. In this case the correction factor was found to 

be 


Pads! ?sio, = 1.57 (27a) 
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Using equations 27 and 27a it is possible to obtain the true 
pressure generated in the gallium arsenide if the pressure 
in the gold layer or the pressure measured by the quartz 
gauge is known. 

The problem of predicting the pressure generated in 
the gallium arsenide by the laser pulse is now reduced to 
predicting the pressure in the gold layer. From equation 11 


it is known that 
2, aay a (28) 
Au Auo p : 


But since the silicon dioxide layer covering the gold layer 
is an absorber of 10.6 um radiation not all of the laser 
energy incident upon the silicon dioxide should be trans- 
mitted to the gold layer. However, from Table III of 
Appendix A, it is seen that the absorption coefficient of 
Sllicon dioxide at 10.6 um is 3.0°10> ena For the silicon 
dioxide films with thicknesses on the order of thousands of 
Angstroms, such as those used in the current investigation, 
the energy fraction transmitted through the silicon dioxide 
Mycrotseapeumt Ul98. So Virtually all of the incident 


energy is transmitted to the gold layer by the silicon dioxide 


layer and it can be stated that 


Bion = Gee (29) 
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where A is the absorptance of the Au-SiO, film system. 


Z 
Equation 28 then becomes 


Pa = (Za, A E,/3t,)*/? (30) 


So if the laser energy density E and the absorptance A 
are known the pressure in the gallium arsenide can be 
calculated using equations 27 and 30. 

2. Derivations of the Absorptance 

Since both the gold and the silicon dioxide layers 
are absorbers of 10.6 um radiation, the problem of finding 
the absorptance of the system is that of finding the 
absorptance of a single absorbing layer (the silicon dioxide) 
on an absorbing substrate (the gold). This problem has been 
solved by Heavens [Ref. 36]. 

The energy that is coupled into the system consists 
of the laser energy that is absorbed by the silicon dioxide 
layer and the laser energy that is transmitted through the 
Silicon dioxide layer and is then absorbed by the gold 
hayern. Lt the absorptance is defined as the ratio of the 
energy coupled into the system by both of these means to 
the incident energy, and the reflectance of the system is 
defined as the ratio of the energy reflected to the energy 


incident, then 


A = 1-R (31) 


Where R is the reflectance. 
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From Heavens, the reflectance of a system of thin 


films such as illustrated in Figure 6 is 


R = cc*/aa* (325) 


where c and a are obtained from the matrix product of the 


matrices of Fresnel coefficients for the system. That is, 


(eee) ja YC (33) 


In general, for an absorbing medium, the elements of the 
matrix of Fresnel coefficients are all complex and are 


given by 


Cea Whee ig xr tis (34) 


eee ae ly 
m Un m m 


where 


eer 
- ie m-1 iil m-1 (35) 
[(mp np )° ‘(my tay 
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E m-1 m m m-1 
ee [(np + Np ) v's (ny + Ny ) ] (36) 
m- m-L 
~ Ber 
“neil 7 we ie (37) 
Y = a d ‘€z:3) 
m-1 ie a m-1 


in which Np and nN; are the real and imaginary parts of the 
complex index of refraction for the material composing the 
layer and d is the thickness of the layer. 

For normal incidence the quantities of equation 34 


are defined by 


pee m-1 cos bane (39) 
met 
qd, @ Sin ey (40) 
r= 2 M1, cos Y hesin yy... ) (41) 
m Sm m-1 m m-1 
sS. = eo Bley cos y Teor sin y. . ) (42) 
m m m-1 om m-1 
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— m-1 e 

ieee = (g,,cos ea h_ sin Nee (43) 
ed 

u, = ¢ (h cos Ym-1 7 Sm522 Yp-7? (44) 
ine ik 

vac cos Y 4 (45) 

eines 
vee -€ sin ¥,_j (46) 


A double suffix notation can be adopted to denote 
the elements of the product matrices. Thus the elements of 


the single layer system under consideration are 


(C1 )€Co) F fPyo * 19% 4 2 * 1542 


f seek BS V + iW 


12 ay a 12 
where 
Pies pe erg, 12 ee 
‘hig © ly" By eile, fae) 
Ea 2 Ep Sys ky, oe) 
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Uj2 = Ay + AYP, - 844, (50) 


and the other terms of the product matrix are obtained in 
the same fashion as equations 47 through 50 by matrix 
multiplication of the individual matrices Coa): From 
equation 32 and from the quantities defined in equations 
35 through 50, the reflectance of the system is 

fete, U5 Pa,,) ($1) 
The absorptance of the system can then be calculated from 
equation 31. 

In general the quantities Np and ny needed to calcu- 
late the terms of equation 51 are wavelength dependent. 
Values for these quantities are tabulated in Table III of 
Appendix A for the wavelength of 10.6 um used in the 


current investigation. 
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IV. EXPERIMENTAL PROCEDURE 


A. MEASUREMENT OF PRESSURE PULSES 

The magnitude of the pressure pulse at the backface 
of the gallium arsenide was monitored using the Sandia 
Quartz Gauge previously described. The output of the 
gauge was displayed aon the screen of a Tektronix 7704 
Oscilloscope and recorded photographically. Triggering for 
this oscilloscope was obtained from a Tektronix T-555 
oscilloscope utilizing its delayed trigger capability. 
The T-555 oscilloscope was, in turn, triggered by the 
current output of the Marx bank of the Carbon Dioxide laser. 
Delays of up to 1.0 microsecond were used to allow for the 
delay between firing the Marx bank and initiation of the 
laser pulse. A schematic of the circuitry used is shown 


ier igure §&. 


B. MEASUREMENT OF ENERGY DELIVERED TO TARGET 

Prior to initiation of the experimental data runs the 
Carbon Dioxide laser was aligned so that its energy output 
at a given applied voltage varied less than ten per cent. 
After the data run was completed the energy output was 
measured using the ballistic thermopile to verify that no 
Significant change in the laser output occurred. The experi- 
mental apparatus was such that continuous monitoring of the 


laser output was not practical. 
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Since the laser beam was required to pass through various 
optical components prior to irradiating the solar cell model 
not all of the energy of the laser pulse reached the target; 
some was absorbed by the various optical components. 

The principal components that absorbed the laser energy 
were the germanium lens used to focus the laser beam and, 
in the case of the vacuum data runs, the window in the test 
chamber through which the laser light was admitted to the 
chamber. 

An experiment was performed in which the energy output 
of the laser was measured with the ballistic thermopile 
both before and after the laser beam passed through the 
germanium lens. The lens was found to transmit about 90% 
of the laser energy. 

A similar experiment was performed on the germanium 
Window used as an entrance window for the vacuum data runs. 
However, aS is described in the Experimental Results section 
of this report, the energy transmitted through the germanium 
window was a function of the number of laser shots. A trans- 
mittance curve was obtained for this window and is shown 
in Figure 9. A correction factor was then applied to the 
laser energy for the value appropriate to the laser shot. 

From the above experiments a correction factor of 0.9 
was applied to the laser energy to account for the absorp- 
tance of the lens and a correction factor corresponding 
to the transmittance of the window for a given laser shot 


was applied to account for the absorptance of the window. 
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V. EXPERIMENTAL RESULTS AND ANALYSIS 


A. ENTRANCE WINDOW DAMAGE 

An unexpected result that arose in the performance of 
mic experiment was the inability to find a material for the 
entrance window to the test chamber that was capable of 
Withstanding the power density of 10/ watts/cm? needed to 
obtain the 1.0 cm spot diameter required on the solar cell 
model. Lack of a window that could withstand more than a 
few shots at this power density severely limited the amount 
and quality of the data obtained in the vacuum data runs. 

Four different materials were tried in an attempt to find 
a suitable material: Sodium Chloride, Potassium Bromide, 
Zinc Selenide, both with and without an anti-reflective 
coating, and anti-reflective coated Germanium. 

Of these materials the Potassium Bromide window showed 
damage in the form of a burn spot on the first shot and the 
Sodium Chloride and the anti-reflective coated Zinc Selenide 
windows showed similar damage after about three shots. 
Closer examination of the Zinc Selenide window showed 
indications of possible internal damage and spallation on 
tess back face. 

Since the anti-reflective coating on the window was 
damaged, it was decided to remove it and determine if the 
uncoated surface was better able to pass the laser beam 
without damage. It was found that this was not the case. 


After three shots, it was observed that further surface 
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damage was occurring and experimentation with this window 
was terminated. 

Finally, an anti-reflective coated germanium window was 
obtained and was found to survive five shots without any 
apparent damage. A vacuum data run was then initiated using 
this window. However, after a few more shots surface damage 
was visible on the window and the run terminated. 

A subsequent data run was made by irradiating a dif- 
ferent portion of the window surface. In this run surface 
damage was observed after five shots. 

Pictures of the damage caused to the Sodium Chloride 
and Germanium windows are included as Figures 10 and ll. 

The primary cause of the damage observed in the window 
materials examined is felt to be that the incident power 
density exceeded the intrinsic threshold for damage of the 
material coupled with damage due to possible inclusions in 
the materials 

The possibility that inclusions in the surface layer of 
the Sodium Chloride and Potassium Bromide windows were 
responsible for the damage to these windows is extremely 
plausible. These windows were polished by the experimenter 
and it is quite likely that small pieces of the polishing 
compound could have been imbedded in the surface of the 
windows. These particles could then have acted as strong 
absorbers for the incident 10.6 um radiation. 

The most probable mechanism to explain the damage 


observed in the Zinc Selenide and Germanium windows is 
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exceeding the intrinsic damage threshold for the anti- 
reflective coatings of these windows. While no information 
for the coating of the Zinc Selenide window was available, 
information given by the manufacturer of the Germanium window 
indicated that the damage threshold for the coating on this 
window was about 100 me gawatts/cm°. Measurement of the burn 
spot on the Germanium window indicated that it had a diameter 
of about 0.5 cm; this corresponds to a power density of 

about 7000 megawatts/cm¢, well in excess of the damage 


threshold of the coating. 


B. PRESSURE PULSE MEASUREMENTS IN VACUUM 

Two data runs were made under a vacuum of approximately 
9° torr. using the germanium entrance window, a laser 
energy of 5.25 joules, and two different sets of thicknesses 
for the silicon dioxide and gold coatings of the solar cell 
model. Tables IV and V of Appendix A list the results for 
each set of film thicknesses. A graph of the natural 
logarithm of the peak pressure measured in the gallium 
arsenide versus shot number of the laser for each set of 
coating thicknesses is displayed in Figures 12 and 14. 

As can be seen on these graphs, the peak pressure in 
each case steadily decreased from a high initial value 
to about 100 bars as the shot number increased. The dura- 
tion of the pressure pulse in both cases remained approxi- 
mately constant at about 400 ns until the last two shots 


when it shortened to around 250 ns. Figures 16 and 17 
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display early and final pressure pulses for the system of 
a 1000 Angstrom layer of silicon dioxide on a 1000 Angstrom 
layer of gold (1000/1000 system). Similar pulses for a 
1000/2000 system are shown in Figures 18 and 19. 
Observations of each gauge after completion of the data 
run revealed that in both cases the layers of silicon 
dioxide and gold had been removed from most of the surface 
of the gallium arsenide substrate. Since the duration of the 
data run for the thicker of the two systems (the 1000/2000 
system) was the shorter (10 shots) it is concluded that it 
required an equal number of shots to remove the coatings 


from the 1000/1000 system. 


C. PRESSURE PULSE MEASUREMENTS AT ATMOSPHERIC PRESSURE 

Three data runs were made at atmospheric pressure using 
the following sets of gold and silicon dioxide layer thick- 
nesses: 

Pee 00/7 2000 

Pee 00 / 5000 

Sey 5000 
In the above sets the thickness of the silicon dioxide layer 
iieenostmomosets tOllowed by the thickness of the gold 
layer in Angstroms. 

In all cases a laser energy output prior to any optical 
components of 5.25 joules was used and a spot diameter of 
1.3 cm was maintained on the solar cell model. The data 
from these runs are tabulated in Tables VI and VII of 


Appendix A. 
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In the 1000/2000 system the peak pressure pulse in the 
gallium arsenide was initially 14.0 bars and decreased to 
a value of 9.0 bars after 13 shots and remained nearly 
constant at this value for the remaining six shots of the 
run. The duration of the pressure pulses measured remained 
nearly constant at 1.0 usec. After six shots, damage to 
the coatings was observed in the form of total removal of 
the silicon dioxide and gold layers at various points on 
the surface of the solar cell model. Figure 20 displays 
a typical pulse for this system. Figure 21 displays a 
graph of the natural logarithm of peak pressure versus shot 
number for this data. 

The second system examined at atmospheric pressure 
consisted of thicker coatings. In this case a 2000 Angstrom 
layer of silicon dioxide was deposited onto a 3000 Angstrom 
layer of gold. The results of the fourteen shots taken at 
this 2000/3000 system did not display the decrease in peak 
pressure observed in the previously described results. For 
this system the peak pressure measured remained nearly 
constant with an average value of 8.0 bars and displayed a 
mean pulse duration of approximately 600 nsec. 

The greater thickness of the silicon dioxide coating 
allowed an observation to be made to obtain an eStimate of 
the number of shots required to at least partially remove 
the silicon dioxide layer. It is estimated that the silicon 
dioxide layer was almost totally removed from the central 


portion of the surface of the solar cell model after only 
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one shot. After a total of six shots, gallium arsenide 
Was Visible at the center of the surface indicating that 
both the silicon dioxide and gold layers were removed. A 
typical pressure pulse measured for this system is shown in 
Figure 22. Figure 24 shows the surface of the solar cell 
model upon completion of fourteen shots. 

A final series of shots was run on a target upon which 
a 5000 Angstrom layer of gold was deposited. No layer of 
Silicon dioxide was deposited on this solar cell model in 
order that the effects of a thick layer of gold upon the 
pressure pulse generated could be examined. An initial peak 
pressure of approximately six bars was measured with a 
duration of approximately 2.0 usec. This magnitude of the 
pressure pulse persisted for four shots and then the pulse 
was lost in the noise displayed on the screen of the 
oscilloscope. The run was terminated at eight shots with 
no damage to the gold surface apparent. Figure 23 shows a 


typical pulse from this series. 


D. COMPUTER MODELING OF PEAK PRESSURE PULSE 

A computer code was developed to predict both the peak 
pressure in vacuum and in atmosphere for a given set of 
Silicon dioxide and gold layer thicknesses. For irradiation 
in vacuum the plasma Benno dei presented by equation 30 
was used and for irradiation at atmospheric pressure the 
laser-supported detonation wave model developed by Pirri and 


expressed by equation 10 was used. 
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In both models the resultant pressure is dependent upon 
the energy density delivered to the target. From equations 
29 and 31 it can be seen that this energy density is at least 
partially dependent upon the reflectance of the target 
irradiated. Consequently, the computer code was designed 
to first calculate the reflectance of the solar cell model 
from equation 32, and from this to obtain the absorptance 
of the system from equation 31. The actual energy density 
delivered to the target was then calculated from equation 29. 
This resultant energy was then applied to equations 10 and 
30 to obtain the expected magnitudes of the pressure pulses 
in vacuum and at atmospheric pressure respectively. For 
the case of the atmospheric measurements the duration of 
the pulse was measured from the photographic records of 
the data and this value was then assigned to the variable 
t in equation 10. 

The code was also designed to examine the effects of 
reducing the thickness of the silicon dioxide and gold 
layers upon the peak pressure expected. Since the obvious 
dependence of the predicted pressure on the thickness of 
the coating layers is through the exponential terms of the 
absorptance, this was accomplished by stripping off the 
Silicon dioxide and gold films in 50-Angstrom increments 
and monitoring the effect this had on the peak pressures 
peedieted. 

The results of running this computer code for the 


thicknesses of silicon dioxide and gold used for the data 
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runs of the current investigation are summarized in Table VIII 
of Appendix A. The computer code itself is given in Appendix F. 
As 1S shown in this table, in all cases where a silicon 
dioxide layer was present the initial reflectance of the 
System was almost 99%. At the point where the entire 
silicon dioxide layer is removed, the reflectance is only 
reduced by about 1%. The value of the reflectance at this 
point corresponds to that of the uncoated gold layer, that 


mmabout 9%. 


E. ANALYSIS OF EXPERIMENTAL RESULTS 
1. Pressure Measurements in Vacuum 
The data depicted in Figures 12 and 14 was fitted 
to a straight line using a linear regression computer code. 
This was done to test the hypothesis that the peak pressure 
measured had a functional dependence upon shot number of 


the form 


eee eee (52) 


where k is some constant dependent upon the initial thicknesses 

of the components of the film system and x is the shot number. 

The linear regression analysis yielded a value for k, the 

peak pressure to be expected after the first shot of the 

laser and a correlation factor that was a measure of the 

degree to which the data fitted the assumed straight line. 
Since the data obtained in the vacuum runs was 


complicated by the fact that the entrance window of the 
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chamber absorbed some portion of the laser energy that 
varied from shot to shot, the measured pressure was 
corrected for the fraction of energy absorbed by the window. 
This was accomplished by dividing the pressure for a given 
shot number by the square root of the transmittance of the 
germanium window for that shot number. The factor qi/2 
was chosen since the energy transmitted through the window 

is TES where ES 1s the energy incident on the window. The 
plasma tamping model used for predicting the pressure in 
vacuum is expressed by equation 30 wherein the energy 

density enters to the one-half power; hence the correction 
Lagtor qi/2 is applied. 

Dividing by this factor has the effect of normalizing 
all the pressures to the values that would have been measured 
had the window not been present. A check to see if the only 
effect the absorption of energy by the window was to lower 
the pressure measured by comparing the value of k calcu- 
mueecdepyertttine a straight line to this "corrected" data 
With that obtained by analysis of the uncorrected data. 
Figures 12 and 14 display a graph of these corrected 
pressures versus shot number. If the values for k were 
equal or nearly so then the dependence of the peak pres- 
sure upon film thickness was independent of the energy 
absorption by the window. 

For the 1000/1000 system the analysis yielded a value 
of k equal to -0.14 before the correction for the window 


was applied; after the correction was applied k was equal 
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mie-02-15,. For the 1000/2000 system the uncorrected value of 
k was -0.15 and after correction was -0.15. Based on the 
good agreement between the values of k before and after 

the correction was applied, the assumption that the energy 
absorption by the window was not dependent upon the thickness 
of the films was considered justified. 

To determine the validity of the hypothesis of the 
functional relationship between shot number and peak pressure 
indicated by equation 52, the correlation factor calculated 
by the linear regression code was examined. In the case 
of the 1000/1000 system a correlation of 0.79 was calculated. 
For the 1000/2000 system the correlation factor was 0.89. 
Since a correlation of 1.0 indicates a perfect straight 
line, the values obtained indicated that the functional 
dependence assumed in equation 52 was possibly justified. 

However, the correlation value of 0.79 calculated 
for the 1000/1000 system seemed to indicate that the 
functional relationship between peak pressure and shot 
number of equation 52 was not entirely correct. Closer 
examination of the data for the 1000/1000 system displayed 
in Figure 11 indicated that the peak pressure measured 
remained relatively constant from shot 10 on. Consequently, 


a functional relationship of the form 


P= P+ Pe (53) 
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where k' was assumed to be thickness dependent was assumed 
for the first nine shots of the 1000/1000 and 1000/2000 
systems. After nine shots it was presupposed that the 

-pressure saturated to P, the average pressure measured for 
shots 10 through 16 of the 100/100 system. Since the data 
for the 1000/2000 system covered only ten shots it was not 
possible to determine a corresponding value of P for that 
System and the value for the 1000/1000 system was also used 
for the 1000/2000 system. 

To verify the relationship suggested by equation 53 
the mean value of the pressures measured for shots ten through 
sixteen of the 1000/1000 system was calculated and found to 
be 215 bar. This value was then subtracted from P, the 
value of peak pressure after correction for the absorptance 
of the germanium window, corresponding to each shot number 
for the results of the vacuum data runs. The results of 
this subtraction are displayed in Tables IX and X of Appendix 
A. A graph of the natural logarithm of the resultant pressure 
P' versus shot number was then made. These graphs are shown 
as Figures 13 and 15. A linear regression analysis was 
performed on these data to test the validity of the 
relationship 

|G pra 


Pi = (ey 2 PD) = ae (54) 


suggested by equation 53 for the first nine shots. 
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For the 1000/1000 system a value of k' of -0.44 
and a correlation of 0.91 were obtained. The corresponding 
results for the 1000/2000 system were k' = -0.22 and a 
correlation of 0.85. Comparing these correlations to those 
meported above for the verification of equation 52 it is 
seen that, in the case of the 1000/1000 system, the correla- 
tion increased significantly. For the other system considered, 
the 1000/2000 system, the correlation did not change 
appreciably. 

Viewing the fact that equation 54 yielded a higher 
correlation for the 1000/1000 system than did equation 52 
and that the results for the 1000/2000 system did not 
appreciably change, it was concluded that the relationship 
between shot number and pressure generated is of the form 
of equation 54. That is, a "saturated" exponential for the 
first nine shots and a Sonar value thereafter. The 
Saturation is felt to occur when the gold and silicon 
dioxide coating are entirely removed. 

An assumption was made that the value of k in 
equation 52 and k' in equation 54 are dependent upon the 
initial thicknesses of the silicon dioxide and gold films 
of the solar cell model. An attempt to at least partially 
justify this assumption was made using the previously 
described computer code. 

If the relationship between peak pressure and f1lm 


thickness exists and is of the form of equation 54, that is, 
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a Saturated exponential dependent upon the initial thicknesses 
of the gold and silicon dioxide films, the dependence will 
appear in the absorptance of the system since that is the 
only place where the film thicknesses appear. If there is 
a relationship between the constant k' in equation 54 and 
the initial thicknesses of the films then the pressure 
predicted by the plasma tamping model should equal that 
measured after the first shot of the laser. 

The initial pressures measured for the film systems 
used (as predicted after the data had been analyzed by the 
linear regression) were: 

1. 1000/1000 system: 1706 bar. 

2. 1000/2000 system: 753 bar. 

The values predicted using equation 30 were: 

1. 1000/1000 system: 556 bar. 

2. 1000/2000 system: 556 bar. 

The predicted value for each caSe is respectively 33% and 
74% of the measured value. Hence, it is concluded that 
there is at least some justification to the assumed rela- 
tionship between film thickness and the proportionality 
constant k' of equation 54. The existence of this rela- 
tionship would then imply that the decrease in peak 
pressure measured for successive shots can be attributed to 
the removal of the gold and silicon dioxide films by suc- 
cessive shots. Unfortunately, the computer model did not 


predict this decrease. 
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2. Pressure Measurements at Atmospheric Pressure 

As has been previously described the measurements 
of pressure pulses generated at atmospheric pressure were 
of small amplitude. These results were as expected since 
previous work by Magee and others [Ref. 29] noted that in 
comparison to pressure pulses measured in vacuum, pressure 
pulses measured in atmosphere on identical targets were 
Significantly smaller. 

In the atmosphere measurements performed in the 
current investigation the pulses measured were generally 
of less than 10 bars magnitude. In only one case, the 
1000/2000 system was there any variation of the pressure 
measured with shot number. The pressures measured on the 
thicker films, namely the 2000/3000 system and the 0/5000 
system, remained almost constant throughout the data run 
with average values of 7.4 and 5.4 bars respectively. 

For the data from the 1000/2000 system a linear 
regression analysis identical to that performed on the data 
from the vacuum cases was performed. However, no corrections 
were necessary for absorption by the entrance window since 
none Was used. The linear regression analysis yielded a 
value of k equal to -0.046 and a correlation to the expo- 
nential relation between pressure and shot number of 0.87. 
This indicates that while an exponential relationship 
between shot number and pressure probably exists, is is a 
rather weak relationship, much weaker than that for the 


Vacuum case. 
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The laser-supported detonation wave model for the 
generation of the pressure pulse leading to equation 10 
was used to attempt to predict the pressure pulse using the 
experimenter-generated computer code previously discussed. 
The results were uniformly dismal. Comparison of the values 
tabulated in Tables VI, VII, and VIII show that the pre- 
dicted pressures were consistently at least two orders of 


magnitude smaller than the experimentally measured pulses. 
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VI. CONCLUSIONS 


A. DEPENDENCE OF PRESSURE UPON SHOT NUMBER AND FILM THICKNESS 
Analysis of the curves of pressure versus shot number 
Shown in Figures 13, 15 and 21 indicates that there exists 
an exponential relationship between the pressure generated 
and the shot number. Since total removal of the silicon 
dioxide and gold films covering the gallium arsenide 
substrate of the solar cell model was noted at some point 
during the data runs it is concluded that there exists some 
relationship between the peak pressure generated and the 
remaining thickness of the thin layers of silicon dioxide 
and gold remaining on the gallium arsenide substrate. The 
fact that the gallium arsenide substrate was entirely 
exposed at some points after only six shots on the 1000/2000 
system at atmospheric pressure and that after this exposure 
was noted the pressure generated remained relatively 
constant seems to validate this conclusion. It is unfortu- 
nate that no similar observations could be made during the 
vacuum cases but the position of the solar cell model in 
the test chamber and the opacity of germanium to visible 
light precluded observation of the damage to the target 


until completion of the data run. 
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B. DETERMINATION OF THE PRESSURE GENERATION MECHANISM 

The primary methods used to attempt to ascertain the 
mechanism whereby the observed pressure pulses were 
generated were examination of the duration of the measured 
pulse and a comparison of the measured pulse magnitude with 
that predicted by whichever of the two models considered 
in the previous sections of this report was applicable to 
the conditions of the particular experiment. - 

The comparison of the initial pressure pulses measured 
in vacuum to those predicted by the plasma tamping model of 
equation 30 yielded reasonable nVengee Nene on the observed 
and predicted pulse amplitudes. However, the consistent 
underestimation of the pulse amplitude by equation 30 indi- 
cates that some other mechanism may possibly be contributing 
Significantly to the generation of the pressure pulses. 

In an attempt to determine what this possible alternative 
mechanism might be an examination of the duration of the 
measured pressure pulse was made. In general, the duration 
of the pulses remained near 400 nanoseconds. This is 
about one and one-half the duration of the 250-nanosecond 
pulse of the Carbon Dioxide laser used to irradiate the 
solar cell model. Yang [Ref. 28] has reported that the 
pressure pulse is expected to have a duration close to that 
of the irradiating laser for pressure pulses generated by 
plasma tamping. This fact gave strong indication that 
possibly some other mechanism than plasma tamping might be 


the dominant mechanism. 
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A further indication to this effect is the requirement 
that power densities on the order of ioe watts/cm- are 
needed for plasma production [Ref. 41]. The power density 
used throughout the current investigation was of the order 


o£ 10” 


watts/cm“, much lower than the power density needed 
for plasma production. This strongly suggested that a 
plasma was not being formed by laser irradiation of the 
target, thereby strongly suggesting that the plasma tamping 
model was not wholly applicable to the conditions of the 
experiment. 

Based on these two observations the plasma tamping model 
was discarded as the chief mechanism for pressure generation 
in vacuum under the conditions of the current investigation. 

The applicability of the laser-supported detonation wave 
model was not considered in the vacuum case since the forma- 
tion of such a detonation wave required pressures far in 
excess of the 1) torr. pressure used in the vacuum measure- 
ments. This model was examined as a possible mechanism for 
the results obtained at atmospheric pressure but again the 
fact that the power density needed for plasma production 
was not met forced its exclusion as a possible mechanism. 
Also, the uniformly poor prediction of the model in the 
computer code would have forced its discard regardless of 
whether or not the conditions for plasma production were met. 

The final mechanism examined as a possible explanation 
of the pressure pulse generation was the thermomechanical 


mechanism described by Ready [Ref. 42]. This mechanism 
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appeared attractive Since it did not require removal of 
material from the target or the formation of a plasma. All 
that is required is rapid enough heating of an absorbing 
medium to cause the thermal expansion which generates the 
pressure wave. The mechanism suggests that the duration of 
the pressure wave generated is not tied to the laser pulse 
length but to the time required for the absorbing medium to 
thermally expand to its greatest extent and than to relax 
to its ambient thickness. In general, this time would 
probably be much longer than the laser pulse. 

The strong possibility that this mechanism is the dominant 
mechanism in the current investigation is supported by two 
factors. The first is the long duration of the pressure 
pulse when compared to that of the laser pulse and the second 
is that the pressure pulse measured for the very thick 0.5000 
System was of much smaller amplitude and greater duration 
than those measured for the thinner 1000/2000 and 1000/1000 
systems. This small amplitude and long duration is attributed 
Bomene ability Of the thick film to absorb and dissipate 
the thermal expansion over a much greater thickness than 
the thin film. This accounts for the long duration of the 
pressure pulse for the thick film since it takes a longer 


time for the thick film to expand than for a thin one. 
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C. DETERMINATION OF POWER THRESHOLD FOR REMOVAL OF 
COVERING FILMS 


In all cases except the 0/5000 system, the layers of 
Silicon dioxide and gold were found to have been rapidly 
removed by the relatively low power density of 107 watts/cm- 
used. throughout the current investigation. In the cases 
where observation was made of the number of shots necessary 
for this complete removal to occur, the atmospheric irradia- 
tions, it was observed that no more than six shots were 
required before gallium arsenide was visible at some point 
on the surface of the solar cell model. 

If the dominant mechanism at work is the thermomechanical 
mechanism described in the previous section, this removal 
of the films of gold and silicon dioxide is attributed to 
Spallation. A portion of the pressure wave generated by 
thermal expansion of the gold layer is reflected at the gold- 
gallium arsenide interface and is of sufficient magnitude 
to exceed the strength of the bond holding the silicon 
dioxide to the gold thereby causing spalling of the 
Silicon dioxide layer. A similar sequence of reflected 
waves could account for the removal of the gold coating. 

It should be noted that this explanation for the removal 
of the covering films was not verified by experimental 
work during the current investigation but is merely offered 


as a possible hypothesis to explain the film removal. 
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D. SUMMARY OF RESULTS 

The main conclusions arrived at in the examination of 
@me results of the current investigation are: 

1. For films of gold less than 5000 Angstroms thick, 
there 1s a relationship between the pressure generated and 


the shot number of the form 


The value of k' is larger in vacuum than for irradiation at 
atmospheric pressure and the value of P in atmosphere is 
probably very much smaller than the corresponding value 

in vacuum. Since removal of the gold and silicon dioxide 
layers was observed after a series of shots, a dependence 

of the constant k' upon the initial thicknesses of the films 
1s indicated. 

2. The principal mechanism for the generation of the 
pressure pulses observed is apparently thermomechanical 
and due to thermal expansion of the gold layer. 

3. The coatings of sideon dioxide used as an anti- 
reflective coating on the solar cell model are highly 
vulnerable to irradiation by a Carbon Dioxide laser even 
at the low power density of ‘ug! watts/cm* used in the 
current investigation. However, this vulnerability may be 
significantly reduced by making the thickness of this layer 


much greater. The small amount of data obtained for the 


67 





irradiation of a 5000 Angstrom thick gold layer indicates 
that the contact grid composed of a gold layer of this 
thickness is able to absorb power densities on the order of 


nO! watts/cm- without damage. 
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VII. RECOMMENDATIONS FOR FURTHER RESEARCH 


Several additional lines of investigation on the effects 
Of laser irradiation on solar cells are indicated by the 
results of the current investigation. The first of these 
areas 1s verification of the results reported herein. 

These conclusions were based on a very small data base and 
hence are more indications of trends rather than fully 
justified arguments. However, before this can be done the 
experimental apparatus requires modification in some way to 
permit data to be taken under vacuum conditions without 
destruction of the entrance window to the test chamber. 
Also, several other sets of film thicknesses should be 
examined since the current investigation shows that the 1000 
and 2000 Angstrom thicknesses of silicon dioxide are too 
thin to have a significant effect on the pressure generating 
mechanism. 

Another possible area of investigation would be an 
attempt to time-resolve the pressure pulse generated, par- 
ticularly to determine its time of onset and to attempt to 
resolve the contributions of each layer coating the gallium 
arsenide to the total pressure pulse measured. 

Finally, since the results of this experiment seem 
to indicate thermomechanical generation of the pressure 
wave, investigation should be conducted to verify that this 


is indeed the mechanism operating. A refinement of the 
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computer model used in the current investigation to accomplish 
this is indicated. Also an attempt to predict the pressure 
pulse as a function of time should be attempted utilizing 


the material presented in Appendix C. 
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APPENDIX A 


TABLES 


TABLE I 


NPS CO. PeAetaser Characteristics 


Energy - - - - - - - 
feesemWidth- - - - - 


Peak Output Power- - 


- 


Wavelength - - - - - 
Beam Shape - - - - - 
Marx Input Voltage - 


Laser Efficiency - - 


a. Only 14.4 joules have been 
of about 21 Megawatt. 


- - -square, 


= “=p to 4s joules® 


=7 510 Anosecondae 


UD Om 2 5 megawatts* 


4.5 cm x 4.5 em 


= 


- - 10.6 um 


e 


66 kilovolt 


-5% minimum 


produced with a peak power 


b. Considers only electron produced peak. 


eeeelicocused, 75 cm from output mirror. 


TABLE 


HADRON Ballistic Thermopile 


Ballistic Time Constant 
Equivalent Noise Level 


Maximum Energy Input - 


Maximum Power Input- 


Sensitivity- - - -- - 


i 


-50 seconds 
microjoules 
1000 joules 


- 10° watts 


- - - - 44,9 microvolts/joule 
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TABLE IV 


Results for Data Run in Vacuum for 


a 1000 A Layer of SiO, on a 1000 A of Au 


2 
Shot Peak Pressure (bar) Pulse Length 
Number Uncorrected Corrected (nsec. ) 
ii = _ age 
2 1440 2100 400 
3 518 1016 300 
4 Zo2 485 400 
5 = ae 2 
6 144 S031 240 
i Daly 452 400 
8 144 3501 400 
r 9 130 LEZ 340 
10 80 166 6410 
ii 97 ZANE: 400 
eZ 97 Ze 400 
WS == a aia 
14 130 272 300 
Us BU 239 280 
16 100 209 320 


igs 
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TABLE V 


Results for Data Run in Vacuum for a 


1000 A Layer of Si0. on a 2000 A Layer of Au 


Z 
Shot Peak Pressure (bar) Pulse Length 
Number Uncorrected Corrected (nsec. ) 
1 -- -- 480 
2 462 764 440 
3 406 oe 440 
4 Sf 661 440 
5 2 | 482 400 
6 Zoo 444 440 
7 = ae 2s 
8 259 540 320 
9 144 501 240 
10 1S 239 280 
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PABLE VI 


Results for Data Run in Atmosphere for 


a 1000 A Layer of Sid, on a 2000 A Layer of Au 


Shot Peak Pressure Pulse Length 
Number (bar) (nsec.) 
il ae ie 
2 14 500 
3 iS 600 
4 ee = 
5 ae — 
6 ae se 
7 14 500 
8 14 500 
9 12 sae 
10 is z= 
1a =e a 
12 -- _ 
3 eZ 600 
14 9 600 
15 9 600 
16 9 500 
ape 7 500 
18 6 500 
19 9 500 


les 





TABLE VII 
Summary of Results for 2000/3000 and 0/5000 


Systems at Atmospheric Pressure 


Total Number Mean Peak Mean Pulse 
System of Shots Pressure Length 
(Bar) (nsec. ) 
2000/3000 14 S210 600 
a 
0/500 7 5.4 2000 


a. Pressure pulses for this system were of such small 
amplitude that it was extremely difficult to separate 
them from system noise. Therefore, the peak pressure 
and pulse duration upon which these mean values are 


based are, at best, estimates of true values. 
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TABLE IX 


Peak Pressure for 1000/1000 System in Vacuum Corrected 


for an Average Saturation Pressure of 215 Bar. 


Shot saturated Peak Peak Pressure Corrected 
Number Pressure (Bar) for Saturation (Bar) 

it 4s ae 

2 2100 1885 
3 1016 801 
4 485 270 
5 a a 

6 301 86 
i 452 Zo 
8 301 86 
9 Dif di ow 
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TABLE X 


Peak Pressure for 1000/2000 System in Vacuum Corrected 


for an Average Saturation Pressure of 215 Bar 


Shot Saturated Peak Peak Pressure Corrected 
Number Pressure (Bar) for Saturation (bar) 

ik = = 

2 764 549 
3 or 576 
4 661 446 
5 482 267 
6 444 229 
7 ae ae 

8 540 a5 
9 Som 86 


ies, 





APPENDIX B 
Calculation of Values of Np» Ny and k for Gallium 


Arsenide, Gold and Silicon Dioxide 


From Ref. 43 the values of the real and imaginary 


parts of the index of refraction are given by 


np = (GH? (B-2) 
where 
: Dae ee : 
8 = 1- (w/o w*)) (B-3) 
| oe (87+ (1-8) 7 (v_/w) 2y}/? (B- 4) 


in which w is the laser angular frequency, W, is the plasma 


frequency given by 


= 
N 


[4ne2n, /m]*/? 


3.5 (pz)!/2-101° (radians/sec) (B-5) 


where 0 is the density and z is the number of free electrons 


per atom. 
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The quantity Va 15) tieweol lisaon @uequency and is 


given by 
v= w*/4n0 (B-6) 
€ p O 
in which - 1s the static conductivity. 
The absorption coefficient, k, is defined by 
k = (2w/c)n, (Bei) 
where c is the speed of light in vacuum 3. ouge” cm/sec). 
The "skin depth" is defined by 
§ =k + = A_/40n (B-8) 
O lt 
For the co, laser, Ne = 10.6 um and w = La Wome radians/ 
ec. 
The transmittance, the ratio of transmitted to 
incident energy, is defined by 
T=e (B-9) 


where x is the thickness of the sample. 
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i Calculation of Np» n; and k for gold. 
For gold z = 1 and a, = 3.7-101” sec™! [Ref. 37]. 


Then 


w = 8.4 (9z)1/2.1015 sec?) 


= 3,7-10/° radians/sec + 


< 
iH 


2 
oe) ae 


piguaig? * zee} 


Then, substituting WH and v. into equations B-3 and B-4, 


it follows 


g=1 - (5 /to® + vey? 


= -1,2°10° 


and 


ZA Z 


+ (1-8)7 (¥,/8)7) 


Beane 


i) 
bo 


Substituting these quantities into equations B-1l and B-2, 
the values for the real and imaginary parts of the index 


Of refraction are 
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Np 14,2 


ny Pes 3S 


From equation B-7, the value of the absorption coefficient 
mor gold is calculated to be k = ieeong cm}, 
gee calculation for Ny and@k tOr Silicon dioxide. 

Reference 53 shows the transmittance of silicon dioxide 


at 10.6am to be 0.82 for a sample 6680 A thick. From 


equation B-9, it is known that 


Tse %* 
Then 
in T = -kx 
and 
k = =in T/x 4 (B-10) 


From the parameters given in Ref. 53, 


k = 3.0°10° cm 


From equation B-/, ny, is found to be 


ny =i 025 
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3. Calculation of Absorption Coefficient for Gallium 
Arsenide. 


From equation B-7, it is known that 
k = (2w/c)ny 


Table III of Appendix A shows that at 10.6 umn 


IL 
Therefore k = 1.0°1077 en. Note that this gives a skin 


depth of about 1.0 meter indicating that gallium arsenide 


essentially transparent to the radiation from a co, laser. 


84 


= 2) Calo 


1S 





APPENDIX C 


Development of Stress (o(x,t)) and Temperature (T(x,t)) 


as Functions of Position, x, and time, t 


An attempt was made to calculate the stress generated 
in the solar cell model as a function of the distance the 
stress has propagated into the model, x, and the time 
elapsed, t. 

Michaels [Ref. 39] has solved for the stress generated 
in a bar as a function of x and t. 

For an elemental slab of a bar of cross-sectional area 


A, the net force acting on the slab is 
= Z 2 
Ado/ox dx = pA dxd u/dt (C-1) 


in one dimension where 


u = longitudinal displacement of slab, 
fece— time, 

x = longitudinal coordinate of bar, 

0 = mass density of bar, 

co = longitudinal stress generated in bar. 


In order for a one-dimensional treatment of the problem 
to be justified the effects of lateral inertia and shear 
must be negligible. The laser spot size used in the 


current investigation irradiated the entire surface of the 
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target so that the effects of lateral inertia and shear 
can be neglected. 


From thermo-elasticity 
o = E(e - aT(x,t)) (e227) 


where 
2 


E = elastic modulus in dynes/cm”, 


0 
a = coefficient of expansion in cm/cm K, 
ie=etemperature of bar, 
eeeecoeca! Aongltudinal strain. 


Since 


€ = du/ox 


Z 


and c° =E/p where c is the stress wave velocity, equation 


C-1 becomes 
g2usax’ - adT(x.t)/3x = 1/c7a-u/at~ . (C-3) 


For a finite bar extending from 0 to L, we have the boundary 
conditions that the stress at each end of the bar 1s zero. 
That 1s 


= holl/ OX a= Cl (out) | = 0 (C-4) 


a | x=0 
x=L 


X 
x 


0 
L 
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or 


3u/ox|,@ 4 = aT(0,t) (C-5) 
du/ ax], = aT(L,t) (C-6) 
The initial conditions for the bar are 


Gieeg 2 2 and du/stiiog = 0 


Recall equation C-3. Taking its finite Fourier cosine transform 


iL 
1/c*3*/at- J US rie] elas Getic l,j) seb’ 
O 


iL 
= J 9 -u/ax? G05 (mic/ i) ce 
O 
iL 
= (6! 4 Olu (Oenuts) 1/6 se OS™ (nll) eacox (C-7) 


fe) 


Integrating the first term on the right-hand side of this 


equation by parts, the result is 


L 
i] a7u/ax? cos (BPX) ax = (-1)P[du/ax],_, - [2u/3x],_, 
O 
L 


néa? nx 
- ee LW (xee eos Spd dx 


0 


Mmivekmae the boundary conditions, equations C-5 and C-6, 


this becomes 
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L 
f 5 
0 


2u/axé cos (==) dx = (ae alti. t) “salt. 4) 


L 
22 
cm f 2,0) cos (“) dx — (C-8) 
0 


Using equations C-5, C-6 and C-8 equation C-7 becomes 


ie- aa (n,t)/até = (-1)"aT(L,t) - aT(0,t) 


where 


and 


- n°n*/L*a_(n,t) - 08 (n,t) (C-9) 
L 
u.(n,t) = J ax,t) cos (ntx/L) dx (e-00p 
O 
L 
@_(n,t) = hs 3T(x,t)/9x cos (ntx/L) dx (eeii 
O 


Bauation C-9 reduces to 


ait (n,t)/at® BF (n?n?c*/L)@_(n,t) = BOG (Cod) EA - T(0,t) 


Now, 


- 0 (n,t)] (C-12) 


define 


ne = Rocce“ /L- and 


P(n,t) = c7a[(-1)"T(L,t) - T(0,t) - 6. (n,t)] 
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and equation C-12 becomes 


Z2— 2 2 = 
3 u/ot * wou, = Ei) COMES) 


The solution to this differential equation is [Ref. 39] 


t 
1/u J P(n,z) sin (w(t-z)) dz 


O 


u.(n,t) 


i 
cLa/nt JOG ane = (0 2) 8 (n,z)] 
O 


Siabol (Gatuge/ Gera) wale (Coe lel) 


By integrating equation C-ll by parts it is found that 


L 
B_(n,t) = T(L,t)(-1)" - T(0,t) + na/L f T(x,t) sin (m x/L) dx 
O 


Substituting this into equation C-14 it is seen that 


ic 1B 
u.(n,t) - -ca f sin(nmc/L(t-z)) J tGx,2) Simi(nwx/ L) dxdz ~(C=15)) 
O O 


byeesine the definition of the finite Fourier transform on 


equation C-15 it is found that 
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T_(0,t) 


cx, t) = Sr - 





3 t L 
2 nwx weet = Z,) : 
cos = ecoulac ees 
L aa L ; TD me Ores z)Sin—— dx 
O 
But 
t L 
u.(0,t) = ca fsin(o) f T(x,z) Sin(0) dxdz = 0 
O O 
So 
m ie L 
x,t) = =2ce 2 cos (AP) f sin(PtE2)) ( 1(x,z)sin (APS) axdz 
L n=l ' 
O O 
(C-16) 


Hence, if T(x,t) is known, the stress in the bar can be 


obtained from equation C-2. That is 
o(x,t) = oc“ [au(x,t)/ax =e Xt) | (C-17) 
The problem of determining T(x,t) for a two-layered 


wall analogous to the silicon dioxide and gold layers of 


the solar cell model used in the current investigation has 


been solved by Griffith and Horton [Ref. 40]. 
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The solution assumed that the area of the wall was 
sufficiently large that the heat is conducted perpendicularly 
to its surface, which is taken as the plane x = 0. The 
wall consists of two components, one of thickness a, extending 
from the plane x = 0 to x = a, while the second component 
extends from x = @ to x = ~, 

The temperature T(x,t) must satisfy the following 


boundary value problem 


g/t = Do-T/ax” £> 0, a >x > 0; D=D, (Gauss 
ce INP ea gar tol ND dD. 

T, (x,0) = 0 (C-19) 

T, (a, t) = T, (a, t) (C-20) 

k, 97, /2x| 2g = K22T2/ 2X1 yay (C=2i) 

(-k, 8T,/8X) | 29 = He Lox Ue scm cee (C-22) 

T, (*, t) =0etore0 = t < (C-23) 


where subscript 1 refers to the top layer l, subscript 2 


refers to the base layers 2 and 
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ct 


Bey ee =, 


S 


distance perpendicular to face of wall, 
time, 

thermal diffusivity = k/S, 

thermal conductivity, 


constant flux of heat perpendicular to the wall 
surface, 


volume specific heat (specific heat times density). 


The solutions to this boundary value problem are found to 


be: 


Dest 2 
T, (x,t) = cet eX /4Dit _ x (1-erf{x/24D,t1) ] 
; 1 


: aa ° - {x- 24/4D 
- : : =)" [s pete eo (nea, /4D,t _ {x-20(n+1)“}/ ) 
Meerenel))}(l-ers (see) 3 
24D t 
+{x-20 PES ler fi (Ca2S). 
| 20D, t 
T. (x,t) = 2HA {Dy : 2h) B¢p| 2° _-{x-a{1- (D,/D; (2n+1)}}°/4D,t 
2% ) ie 5 yf T 


D x-a{1-{D Paes mets) | 
- [x-a{1-q 2? ey ser. (eae) 
1 


24Dqt 


The surface temperature Wa Oe) isscivien by setting x = 0 


in equation C-24. Then 


OZ 





Die t 00 Dt 
| = U 4 Si somer tee Bae +] ¢7/tD 
a) = H/k, [2V—— Dp ea IS 


2 a(nt1) (1-er£[2B))) }) (C-25) 
Ni yt 


As can be seen, by applying equations C-24 and C-25 
to equation C-17 it is, in principle, possible to predict 
the stress in the solar cell model as a function of both 
position and time. A computer code to accomplish this is 


available from Dr. Tom Magee of Stanford Research Institute. 
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APPENDIX D 


NPS co, Laser Operating Procedure 


The following recommended procedure includes changes 
made to the operating procedure given by Strickland 
[Ref. 30] and is included as a reference for future users 
of the laser system. The procedure is as follows: 

1. Open main valves at gas tanks and set flow rate 


at following values: 


He = 4.0 
co. = 5.0 
N, = 4,0 


Allow a flow of four minutes before firing the laser. 

2. Set trigger generator to warm-up position and allow 
a two-minute warm-up period. Then set trigger generator 
to single shot position. 

Bem turn interlock switch on Circuit 8 to the "ON" 
position. This activates the red warning lights at entrances 
to the laboratory and puts power to the DC power supply of 
the laser. 

4. Ensure that the north door to the laboratory is 
locked to prevent unwary individuals from entering by this 
means. 

S$. Turn on DC power supply utilizing the switch located 
on front cover. This will light both the amber and green 
lights on the power supply. Allow a two-minute warm-up 


period. 
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6. Put on safety goggles with plexiglass window. 

7. Set Marx generator spark gap pressure corresponding 
to desired firing voltage following posted chart. 

8. Set desired Marx generator discharge voltage on 
variac at the DC power supply. 

9. Set Marx generator charging voltage on preset 
voltage meter at the DC power supply. The voltage set 
Should be one-third of the Marx discharge voltage. 

10. Initiate charging of the system by pressing the 
remote charge switch mounted next to the trigger generator 
or by pressing the ''CHARGE" switch at the DC power supply. 
This will extinguish the green light and activate the red 
light on the power supply. 

11. When the preset voltage level is reached, the audible 
alarm will be activated and the red light will be extinguished. 
Firing of the laser can now be accomplished by pushing the 
button on the trigger generator. If the laser is not fired 
within five seconds after completion of charging the auto- 
matic dump operates and the system is discharged to ground. 

12. If necessary, the system can be dumped at any time 
after initiation of charging either with the remote dump 
switch at the trigger generator or the local dump switch 
at the DC power supply. 

13. Further firing at the preset voltage can be 
accomplished by repeating the procedure beginning at step 8 


or at another voltage beginning at step 5. 
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APPENDIX E 


NPS CO, Laser Alignment Procedure 


The following alignment procedure is suggested for 
future users of the co, Laser System. It should be imple- 
mented whenever the output of the laser is found to have 
significantly changed or when the Sodium Chloride (NaCl) 
Brewster Windows are installed after polishing. The 
alignment procedure requires, in addition to the CO, 
laser itself, the following equipment: 

1. Helium-Neon (HeNe) CW laser with laser aligner 
attachment. 

Zi CO. laser alignment blocks (2) 

3. HADRON Ballistic Thermopile 

eee id bamicrovoltmeter 

5. Copper Mirror (f = 1.0m). 

The alignment procedure assumes that the NaCl windows are 
not installed at initiation of the procedure. It is broken 
into two sections: rough alignment and fine tuning. 
Referring to Figure 25, the procedure is as follows: 

A. Rough Alignment 

1. Remove Germanium exit mirror and its mount. 

2. Install laser alignment blocks. 

3. Install NaCl windows. 


4. Mount HeNe laser about 1.0 meter in front of output 


end of CO. laser. 
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5. Adjust beam from HeNe laser so that it passes 
through the holes in center of laser alignment blocks. 
6. Place a sheet of Mylar between output end of CO, 
laser and HeNe laser. 
7. Adjust Copper Mirror (f£ = 3.0 m) at back end of 
CO, laser resonant cavity until the beam from the HeNe 
laser reflected by this mirror passes back through holes 
in alignment blocks onto the Mylar sheet. 
8. Adjust the Copper Mirror until the reflected beam 
is coincident with the beam from the HeNe laser incident 
on the Mylar. 
9. Reinstall Germanium output mirror and mount. 

10. Adjust germanium mirror until the HeNe beam reflected 
from this window onto the Mylar sheet is coincident with 
the incident beam. 

il. The cO., laser is now rough aligned. Remove the 
laser alignment blocks. This will require removal of the 
NaCl windows. When removing these windows do not change 
their orientation. 

12. Reinstall NaCl windows. 

13. Remove HeNe laser and Mylar sheet. 

14. Place a sheet of laser "footprint" paper about 
20 cm in front of the germanium output window. 

iSeeeeibieht off the CO, laser system. 

16. Starting at about 30 kV discharge voltage from Marx 


generator fire the laser. If no beam footprint is observed 
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on the footprint paper, increase the discharge voltage in 
increments of 3 kV until a footprint is seen. Do not 
exceed 45 kV at this stage of alignment. When footprint of 
laser beam is seen proceed to step 17. 

17, Remove footprint paper and begin fine tuning the 
aeser. 

B. Fine Tuning 

18. Install 1.0 meter focal length copper mirror about 
1.0 meter from output window of laser. 

19. Connect input lead of voltmeter to calibration 
terminals of ballistic thermopile. | 

20. Mount ballistic thermopile about 1.0 m from copper 
mirror of step 18. 

21. Adjust copper mirror so that beam from CO. laser 
enters INPUT cone of ballistic thermopile. 

22. Using the voltmeter to monitor the total laser 
energy output adjust the mirror at back end of laser 
resonant cavity until laser energy output is at a maximum. 
A wait of about two minutes between laser firings will 
yield best results. 

23. Repeat step 22 with germanium output mirror of co, 
laser. This completes the alignment procedure. A calibra- 
tion curve for laser output as a function of discharge voltage 


may now be obtained. 
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APPENCIX F 
COMPUTER CCDE TO CALCULATE PEAK FRESSURE PULSES 
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